Bootstrapping phylogenies inferred from rearrangement data

ABSTRACT: Large-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models. We present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches. Our method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.     

Contribution of sucrose synthase, ADP-glucose pyrophosphorylase and starch synthase to starch synthesis in developing pea seeds

Using genetic variability existing amongst nine pea genotypes (Pisum sativum L.), the biochemical basis of sink strength in developing pea seeds was investigated. Sink strength was considered to be reflected by the rate of starch synthesis (RSS) in the embryo, and sink activity in the seed was reflected by the relative rate of starch synthesis (RRSS). These rates were compared to the activities of three enzymes of the starch biosynthetic pathway [sucrose synthase (Sus), ADP-glucose pyrophosphorylase and starch synthase] at three developmental stages during seed filling (25, 50 and 75% of the dry seed weight). Complete sets of data collected during seed filling for the nine genotypes showed that, for all enzyme activities (expressed on a protein basis), only Sus in the embryo and seed coat was linearly and significantly correlated to RRSS. The contribution of the three enzyme activities to the variability in RSS and RRSS was evaluated by multiple regression analysis for the first two developmental stages. Only Sus activity in the embryo could explain, at least in part, the significant variability observed for both the RSS and the RRSS at each developmental stage. We conclude that Sus activity is a reliable marker of sink activity in developing pea seeds.     

Sea urchin Hox genes: insights into the ancestral Hox cluster

We describe the Hox cluster in the radially symmetric sea urchin and compare our findings to what is known from clusters in bilaterally symmetric animals. Several Hox genes from the direct-developing sea urchin Heliocidaris erythrogramma are described. CHEF gel analysis shows that the Hox genes are clustered on a < or = 300 kilobase (kb) fragment of DNA, and only a single cluster is present, as in lower chordates and other nonvertebrate metazoans. Phylogenetic analyses of sea urchin, amphioxus, Drosophila, and selected vertebrate Hox genes confirm that the H. erythrogramma genes, and others previously cloned from other sea urchins, belong to anterior, central, and posterior groups. Despite their radial body plan and lack of cephalization, echinoderms retain at least one of the anterior group Hox genes, an orthologue of Hox3. The structure of the echinoderm Hox cluster suggests that the ancestral deuterostome had a Hox cluster more similar to the current chordate cluster than was expected Sea urchins have at least three Abd-B type genes, suggesting that Abd-B expansion began before the radiation of deuterostomes.      

Smoking cessation and bronchial epithelial remodelling in COPD: a cross-sectional study

Background

Chronic Obstructive Pulmonary Disease (COPD) is associated with bronchial epithelial changes, including squamous cell metaplasia and goblet cell hyperplasia. These features are partially attributed to activation of the epidermal growth factor receptor (EGFR). Whereas smoking cessation reduces respiratory symptoms and lung function decline in COPD, inflammation persists. We determined epithelial proliferation and composition in bronchial biopsies from current and ex-smokers with COPD, and its relation to duration of smoking cessation.

Methods

114 COPD patients were studied cross-sectionally: 99 males/15 females, age 62 ± 8 years, median 42 pack-years, no corticosteroids, current (n = 72) or ex-smokers (n = 42, median cessation duration 3.5 years), postbronchodilator FEV₁ 63 ± 9% predicted. Squamous cell metaplasia (%), goblet cell (PAS/Alcian Blue⁺) area (%), proliferating (Ki-67⁺) cell numbers (/mm basement membrane), and EGFR expression (%) were measured in intact epithelium of bronchial biopsies.

Results

Ex-smokers with COPD had significantly less epithelial squamous cell metaplasia, proliferating cell numbers, and a trend towards reduced goblet cell area than current smokers with COPD (p = 0.025, p = 0.001, p = 0.081, respectively), but no significant difference in EGFR expression. Epithelial features were not different between short-term quitters (<3.5 years) and current smokers. Long-term quitters (≥3.5 years) had less goblet cell area than both current smokers and short-term quitters (medians: 7.9% vs. 14.4%, p = 0.005; 7.9% vs. 13.5%, p = 0.008; respectively), and less proliferating cell numbers than current smokers (2.8% vs. 18.6%, p < 0.001).

Conclusion

Ex-smokers with COPD had less bronchial epithelial remodelling than current smokers, which was only observed after long-term smoking cessation (>3.5 years).

Trial registration

NCT00158847     

In planta production of a candidate vaccine against bovine papillomavirus type 1

Bovine papillomavirus type 1 (BPV-1) is an economically important virus that induces tumourigenic pathologies in horses and cows. Given that the BPV-1 L1 major coat protein can self-assemble into highly immunogenic higher-order structures, we transiently expressed it in Nicotiana benthamiana as a prelude to producing a candidate vaccine. It was found that plant codon optimization of L1 gave higher levels of expression than its non-optimized counterpart. Following protein extraction, we obtained high yields (183 mg/kg fresh weight leaf tissue) of relatively pure L1, which had self-assembled into virus-like particles (VLPs). We found that these VLPs elicited a highly specific and strong immune response, and therefore they may have utility as a potential vaccine. This is the first report demonstrating the viable production of a candidate BPV vaccine protein in plants.     

Chemopreventive and renal protective effects for docosahexaenoic acid (DHA): implications of CRP and lipid peroxides

Background

The fish oil-derived ω-3 fatty acids, like docosahexanoic (DHA), claim a plethora of health benefits. We currently evaluated the antitumor effects of DHA, alone or in combination with cisplatin (CP) in the EAC solid tumor mice model, and monitored concomitant changes in serum levels of C-reactive protein (CRP), lipid peroxidation (measured as malondialdehyde; MDA) and leukocytic count (LC). Further, we verified the capacity of DHA to ameliorate the lethal, CP-induced nephrotoxicity in rats and the molecular mechanisms involved therein.

Results

EAC-bearing mice exhibited markedly elevated LC (2-fold), CRP (11-fold) and MDA levels (2.7-fold). DHA (125, 250 mg/kg) elicited significant, dose-dependent reductions in tumor size (38%, 79%; respectively), as well as in LC, CRP and MDA levels. These effects for CP were appreciably lower than those of DHA (250 mg/kg). Interestingly, DHA (125 mg/kg) markedly enhanced the chemopreventive effects of CP and boosted its ability to reduce serum CRP and MDA levels. Correlation studies revealed a high degree of positive association between tumor growth and each of CRP (r = 0.85) and leukocytosis (r = 0.89), thus attesting to a diagnostic/prognostic role for CRP. On the other hand, a single CP dose (10 mg/kg) induced nephrotoxicity in rats that was evidenced by proteinuria, deterioration of glomerular filtration rate (GFR, -4-fold), a rise in serum creatinine/urea levels (2–5-fold) after 4 days, and globally-induced animal fatalities after 7 days. Kidney-homogenates from CP-treated rats displayed significantly elevated MDA- and TNF-α-, but reduced GSH-, levels. Rats treated with DHA (250 mg/kg, but not 125 mg/kg) survived the lethal effects of CP, and showed a significant recovery of GFR; while their homogenates had markedly-reduced MDA- and TNF-α-, but -increased GSH-levels. Significant association was detected between creatinine level and those of MDA (r = 0.81), TNF-α ) r = 0.92) and GSH (r = -0.82); implying causal relationships.

Conclusion

DHA elicited prominent chemopreventive effects on its own, and appreciably augmented those of CP as well. The extent of tumor progression in various mouse groups was highly reflected by CRP levels (thus implying a diagnostic/prognostic role for CRP). Further, this study is the first to reveal that DHA can obliterate the lethal CP-induced nephrotoxicity and renal tissue injury. At the molecular level, DHA appears to act by reducing leukocytosis, systemic inflammation, and oxidative stress.